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Abstract 

In this paper, we discuss the minimal number rj of observables Qi, . . . , Qr,, 
where expectation values at some time instants t\, . . . ,tr determine the 
trajectory of a d- level quantum system ("qudit") governed by the Gaus- 
sian semigroup 
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We assume that the macroscopic information about the system in ques- 
tion is given by the mean values £j{Qi) — Tr {Qip{tj)) of n selfadjoint 
operators Qi,...,Qn at some time instants ti < t2 < . ■ . < tr, where 
n < — 1 and r < deg/i(A,L). Here /i(A,L) stands for the minimal 
polynomial of the generator 

U = -\[H,[H,p]] 

of the Gaussian flow 



1 Introduction 

The ability to create, manipulate and characterize quantum states is becoming 
increasingly important in physical research with implications for other areas 
of science, such as: quantum information theory, quantum communication and 
computing. According to one of the basic assumptions of quantum mechanics, 
the achievable information about the state of a physical system is encoded in the 
density matrix p, which allows one to evaluate all possible expectation values of 
observables trough the formula 

(g) - Tr(pQ), (11) 



1 



where Q is a self-adjoint operator representing a particular physical quantity. 

Thus, in order to have full information about a given quantum system we 
need to know its density matrix p. In particular, a very useful tool in this regard 
is quantum state tomography (QST) which provides means for the complete 
reconstruction of the density matrix for a qudit (or a set of qubits) . The general 
procedure relies on the ability to reproduce a large number of identical states 
and to perform a series of measurements on complementary aspects of the state 
within an ensemble. 

Suppose that we can prepare a quantum system repeatedly in the same state 
and make a series of experiments such that we can measure the expectation 
values 

StiQ) = Tr(Qp(t)) (12) 

of some observables Qi, . . . , Qn at different time moments ti < t2 < ■ ■ ■ < tr- 
The fundamental question arises: 

Can we find the average value of any desired operator from the set of measured 
outcomes of a given set Qi, . . . , Q„ 

^ti {Qn) ■ • • ^tr {Qn) 

where < ti < . . . < tr < T , for an interval [0, T] ? 

Among the existing tomographic techniques for quantum systems, the so- 
called homodyne tomography has received much attention in the literature 
PmiSl. In the phase -space formulation of quantum mechanics there is a one- 
to-one relation between a quantum state and the so-called Wigner function. 
Its marginals are accessible experimentally, and an inverse (Radon) transforma- 
tion allows one to reconstruct the phase-space distribution associated with the 
unknown quantum state. 

The question of how to reconstruct states of d-level systems (qudits) is also 
natural. In this case various methods have been proposed to determine p ^ 
|S1 If the problem under consideration is static, then the state of a d-level 
open quantum system (a qudit) can be uniquely determined only if n — d^ — 
1 expectation values of linearly independent observables are at our disposal. 
However, if we assume that we know the dynamics of our system, i.e. we know 
the generator of the time evolution, then we can use the stroboscopic approach 
based on (|13|) . In general, the term "tomography" will be used collectively to 
denote any kind of state-reconstruction method. 

With reference to the terminology used in the system theory, we assume the 
following definition: 

Definition 1 A d-level open quantum system S is said to be {Qi, ■ ■ ■ ,Qn)- 
reconstructible on an interval [0,T], if there exists at least one set of time 
instants < ti < . . . < tr < T such that the state trajectory can be uniquely 



(13) 
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determined by the correspondence 

[0, T] 3 t, ^ £t^ - Tr (p(t,)gO (14) 
for I = 1, . . . ,n, j = 1, . . . ,r. 

The above definition is equivalent to the foUowing one 

Definition 2 A d-level open quantum system S is said to be {Qi, ■ ■ ■ ,Qn)- 
reconstructible on an interval [0,r], if for every two trajectories with distinct 
initial states there exists at least one i e [0,T] and at least one operator Qk G 
{Qi, . . . , Qn} such that 

Tr(QfePi(£)) ^ Tr(Q,p2(i)). (15) 

Remark 1 In the above definitions we assume that the time evolution of the 
system is given in terms of a completely positive semigroup of operators. Ar- 
guments in favour of completely positive semigroups as the foundation of non- 
Hamiltonian dynamics as well as the discussion of properties of such semigroups 
can be found in papers of Kraus [7] , Lindblad [S] , and Gorini et al. |S| . 

In particular, in Lindblad's paper (S] the general form of the generator of 
an arbitrary completely positive semigroup was derived. A linear operator L 
on a set of linear operators BiTi.) := M{'C^), where Ti. ~ is a d-dimensional 
Hilbert space and M denotes the set of matricies with complex entries, proves 
to be the generator of a completely positive semigroup if and only if it can be 
represented in the form 

where Vj G B{T-C) for j = 1, . . . , k, and _ff is a self-adjoint operator also belonging 
to B{n) (cf. also 10 ). 

Remark 2 It is important that for the number r of time instants ti, . . . ,tr we 
do not formulate any restriction (except that it is finite). 

Remark 3 The question of obvious physical interest is to find the minimal 
number of observables Qi, . . . , Q,, for which the rf-level quantum system S with 
the generator L can be (Qi, . . . , (5,,)-reconstructible. It can be shown that if 
the time evolution of the system S is described by the master equation, 

|p(i)=MO, (17) 
then there exists a set of observables Qi, ■ ■ ■ ,Qrj, where 

ri = max (dimKer (All - L)| (18) 
Ae(T(L) I- J 

such that the system S is {Qi, ■ ■ ■ , (3^)-reconstructible. Moreover, if we have 
another set of observables Qi^ . . . ,Qfi with this property, then fj > rj. The 
number rj given by H18|l we will call the index of cyclicity of the system S. 



3 



2 Polynomial Representation of Flows 



The main idea of the stroboscopic approach to quantum tomography is based 
on the polynomial representation of the flow defined by the general master 
equation. Namely, we have 

m — 1 

- exp(Li) = Y.ak{t)]L\ (21) 

A:=0 

where by Cauchy's theorem 

In the above expression dD is any simple closed contour enclosing the spectrum 
of the operator L in the complex plane and 

m— 1 
fe=0 

denotes the minimal polynomial of the generator L. It is interesting that there 
are ways to compute the functions ak{t) in H21|l without summing the expo- 
nential series or without knowing the Jordan canonical form of L. Namely, 
differentiating (|21|l with respect to t and using the minimal polynomial of L 
one finds that the functions ak(t) for k = 0, . . . , m — 1 satisfy the system of 
equations 



dt 
dai{t) 
dt 

da^_i{t) 
dt 



= dQarn-l{t) , 

= aoit) + diam-iit) , (24) 



with initial conditions a^*''(0) = Sik- It can be shown that functions ak{t) are 
mutually linearly independent, therefore for a given T there exists at least one 
set of time instants ti, . . . ,tm {m = deg/x(A,L*)) such that < < . . . < tm < 
T and det [afc(ij )] ^ 0. 

Taking into account the above conditions one finds that the state p(0) can be 
determined uniquely (and the trajectory <^{t)p{0) can be reconstructed) if and 
only if operators of the form {L,*)^Qi for j — 1, . . . , n and fc = 0, 1, . . . , m — 1 
span the space B{H). 

If the dynamical semigroup is completely positive, then the general form of 
the generator L is given by (|16|l . In this case the criterion for reconstructibility 
of a d- level quantum system can be formulated using the operators H and Vj. 
In particular, if we consider an isolated quantum system characterized by a 
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Hamiltonian Hq and Vj = for j = then the minimal number of 

observables Qi, . . . , Qn for which the system is {Qi, . . . , Q,,)-reconstructible is 
given by 

77 = nj+nl^ ^nj, 

where rii = dimKer(Ai/ — Hq) for all Xi E (t{Ho), i = l,...,r (for details 
cf. 0). 

3 The Minimal Number of Observables for Qu- 
dits Governed by Gaussian Semigroups 

Let us assume that the time evolution of a d-level quantum system S is described 
by the generator L given by 

hp = ^{[Hp,H] + [H,pH]}= -^[H,[H,p]] (31) 

that is, the semigroup ^{t) has the form (cf. e.g. pT] ) 

cso 

^ f dse-"'/(2t)e-'ff«pe'^^ (32) 

\/2'Kt J 

—00 

The symbol H in H31|) and (|32|l denotes a self-adjoint operator with the spectrum 

aiH) = {Ai,...,A„}. (33) 

In the sequel rii stands for the multiplicity of the eigenvalue Xi for i — 1, . . . , m. 
One can assume that the elements of the spectrum of H are numbered in such 
a way that the inequalities Ai < A2 < . . . < Am are fulfilled. The following 
theorem holds: 

Theorem 1 The index of cyclicity of the Gaussian semigroup with a generator 
L given by is expressed by the formula 

r] = max{K,7i, . . . ,7^} , (34) 

where r — (jn — l)/2 if m is odd or r — (jn — 2)/2 if m is even, and 

K nl+nl + ... + nl^, (35) 

7n—k 

7fc 2^ Hi Ui+k ■ (36) 

Proof. In order to determine the value of 77 for the generator IL defined by 131|l 
we must find the number of nontrivial invariant factors of the operator L. Let 
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us observe that if a(H) = {Ai, . . . , A,„} then the spectrum of the operator L is 
given by 

cr(L) = ji^y e M ; i^ij = (Ai - Aj)^ , z, j 1, . . . , to | . (37) 

The above statement foUows from the fact that the operator L can also be 
represented as 

L = (g>l+l(g> H"^ ~2H (g) H , (38) 

where 1 denotes the identity in the space B{Ti). Since H is self-adjoint therefore 
the algebraic multiplicity of A^, i.e. the multiplicity of A^ as the root of the 
characteristic polynomial of H, is equal to the geometric multiplicity of A^, 
rii = dimKer (A^l — H) . Of course, we have ni + . . . + rim = dimTi. 

From (|38|) we can see that the multiplicities of the eigenvalues of the operator 
L are not determined uniquely by the multiplicities of Ai € o'(-ff). But if we 
assume that Ai < . . . < Am and Afe — {k — l)c + Ai, where fc = 1, . . . , m, and 
c = const > 0, then the multiplicities of all eigenvalues of L are given by 

7|,_j| = dimKer [(A, - Aj)2]l-L] (39) 

for i ^ j and 

dimKer(L) ^ nl + . . . + n^^ = k (310) 

when i — j. Now, as we know, the minimal number of observables Qi, . . . , Qrj 
for which the qudit S can be (Qi, . . . , Q,,)-reconstructible is given by H18|l . so 
in our case 

T] = max (dimKer[(Aj-Aj)2]l-L]|, (311) 

i.j — l.....m, L J 

where A^ € o'(i?). Using the above formulae and the inequality jk < k for k > r, 
where r is given by (to — l)/2 if to is odd and (to — 2)/2 if to is even, we can 
observe that also without the assumption A^ = (fc — l)c + Ai one obtains 

77 = max{K,7i, . . . ,7^} . (312) 

This completes the proof. 

According to Theorem ^ if the quantum system governed by a Gaussian semi- 
group is (Qi, . . . , (5„)-reconstructible then the number n of observables must 
satisfy the inequality n > rj. In this case there exists a set of time instants 
ti < t2 < ■ ■ ■ < {m = deg/i(A,L)) such that the knowledge of the expec- 
tation values £j{Qi) ~ Tr {p{tj)Qi) for i — 1, . . . , n and j — 1, . . . , to uniquely 
determines the trajectory of the system. 

The natural question arises: what are the criteria governing the choice of 
time instants ti, . . . , trrJ The following theorem holds: 

Theorem 2 Let us assume that < ti < t2 < ■ ■ ■ < t.„i < T . Suppose that the 
mutual distribution of time instants ti, . . . ,t„i is fixed, i.e. a set of nonnegative 
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numbers ci < . . . < c„,. is given and tj := Cjt for j — 1, . . . ,m, and t G . 
Then for T > the set 

t(T) := |(ti,...,t„0 : t,=Cjt, 0<i<— I 

contains almost all sequences of time instants ti, . . . ,tm, i-c. all of them except 
a finite number. 



Proof. As one can check, the expectation values £tj (Qi) and the operators 
(L*)'^Qii are related by the equality 

m — 1 

St.iQi) = ^a4c,t)((L*)'=Q„po), (313) 

k=0 

where we assume that tj = cjt and the bracket (•, •) denotes the Hilbert-Schmidt 
product in B{Ti). One can determine po from H313|) for all those values t G M_|_ 
for which the determinant a{t) is different from zero, i.e. 

a{t) := det[ak{cjt)] ^ 0. (314) 

One can prove that the range of the parameter t e M+ for which a{t) — 
consists only of isolated points on the semiaxis M+, i.e. does not possess any 
accumulation points on M_|_. To this end let us note that since the functions 
t —f ak{t) for k = 0, l,...,m — 1, are analytic on R, the determinant a{t) 
defined by (|314|l is also an analytic function of f G M. If a{t) can be proved to 
be nonvanishing identically on M, then, making use of its analyticity, we shall 
be in position to conclude that the values of t, for which a{t) = 0, are isolated 
points on the axis R. 

It is easy to check that for k = m{m — l)/2 



d'^ait) 



dt^" 



n (c^-c^)- (315) 

l<j<i<m 



According to the assumption ci < C2 < . . . < c™, we have a'^'^^Q) 7^ if 
k = m(m— 1)/2. This means that the analytic function t — > a{t) does not vanish 
identically on M and the set of values of t for which a{t) — cannot contain 
accumulation points. In other words, if we limit ourselves to an arbitrary finite 
interval [0,T], then a{t) can vanish only on a finite number of points belonging 
to [0,T]. This completes the proof. 
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